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摘要 : 在 地 球 同 步 轨 道 特征 的 基础 上 ,介绍 了 地 球 同 步 轨 道 空间 碎片 的 运动 特征 。 针 对 
地 球 同步 轨道 暗 弱 空间 碎片 进行 观测 的 搜索 方案 ， 考 虑 目标 的 预报 误差 以 及 视 场 的 大 小 ， 用 
数学 方法 分 别 计算 了 同心 圆 搜 索 与 仆 楼 梯 搜 索 的 成 功 概率 ， 并 由 计算 结果 分 析 了 两 者 的 特 
征 ， 认 为 在 一 般 情 况 下 ， 同 心 圆 搜索 方案 具有 更 高 的 搜索 效率 。 
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空间 碎片 又 称 为 轨道 碎片 ， 是 指 宇 宙 空 间 中 除了 正常 工作 的 飞行 器 外 所 有 的 人 造物 体 ， 大 到 报废 
的 卫星 整体 ， 小 到 发 动机 点 火 产生 的 粉 未 ， 都 属于 空间 碎片 的 范畴 5" 。 从 1957 年 10 月 4 日 前 苏联 第 
1 颗 人 造 卫 星人 轨 ， 人 类 进入 了 航天 时 代 。 在 过 去 的 半 个 多 世纪 里 ， 人 类 已 经 进行 了 4 000 多 次 航天 
发 射 活动 ， 将 6 000 多 个 航天 器 送 入 了 宇宙 空间 2 。 目 前 ， 具 有 独立 发 射 能 力 的 国家 已 超过 10 个 ， 
拥有 航天 资产 的 国家 超过 50 个 ， 在 用 的 航天 器 数量 近 千 个 。 航 天 活动 带 来 诸多 好 处 的 同时 也 带 来 了 
很 多 负面 影响 ， 其 中 严重 的 负面 影响 是 空间 碎片 的 不 断 增 加 。 据 统计 ， 地 球 轨道 上 可 编目 的 空间 物体 
数量 (尺度 在 10 cm 以 上 ， 绝 大 部 分 是 碎片 ) 从 1981 年 的 5 000 个 增长 到 了 2016 年 4 月 的 17 000 多 
个 2( 图 1) ， 从 图 1 可 以 看 出 ， 这 些 物体 中 解体 产生 的 碎片 (红色 线条 ) 占 了 很 大 一 部 分 。 
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图 1 空间 碎片 数量 变化 


Fig.1 Variation of space debris quantities 
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空间 碎片 绝 大 部 分 在 地 球 低 轨道 (Low Earth Orbit，LEO) 上 ， 另 外 有 相当 一 部 分 在 地 球 同 步 轨道 
(Geostationary Earth Orbit,， GEO) 上 (如 图 2)。 虽 然 地 球 同步 轨道 碎片 的 密度 远 小 于 低 轨 的 密度 ,但 
地 球 同步 轨道 没有 类 似 大 气 阻 尼 这 样 的 碎片 清除 机 制 ” ， 导 致 地 球 同 步 轨道 的 碎片 在 轨 时 间 长 ， 数 
目 增长 快 ， 对 人 类 航天 活动 具有 潜在 的 危险 。 地 球 同步 轨道 是 一 种 非常 珍贵 的 资源 ， 该 区 域 的 空间 碎 
片 引起 了 航天 国家 的 高 度 关注 。 
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Fig.2 Distribution of space debris density 
(红色 : 直径 大 于 1 mm 的 碎片 ; 绿色 : 直径 大 于 1 cm 的 碎片 ; 蓝 色 : 直径 大 于 10 cm 的 碎片 ) 


地 球 同步 轨道 高 度 为 35 785 km ， 由 于 雷达 的 探测 信号 与 一 成 正比 ”， 地 基 雷 达 一 般 用 于 地 球 低 
轨道 上 空间 目标 的 测量 ， 很 难 实现 对 地 球 同 步 轨 道 空间 碎片 的 探测 。 因 此 ， 从 20 世纪 90 年 代 开 始 ， 
世界 上 主要 的 航天 大 国 先后 将 地 面 光 学 设备 应 用 于 地 球 同步 轨道 空间 碎片 的 观测 研究 ， 地 球 同 步 轨 道 
空间 碎片 的 光学 观测 主要 包括 目标 检测 与 识别 、 目 标的 精密 定位 与 轨道 确定 以 及 目标 的 编目 处 理 等 。 
近年 来 ， 将 空间 目标 视 为 非 质点 ， 通 过 测 光 观测 以 及 多 波段 观测 ， 研 究 空间 目标 的 相关 物理 特性 成 为 
主要 研究 内 容 之 一 后。 为 了 对 目标 成 功 进行 测 光 ， 需 要 在 望远镜 导 星 镜 视 场 内 成 功 地 搜索 并 捕获 目 
标 ， 因 此 高 效 的 搜索 策略 非常 重要 。 


1 地 球 同步 轨道 空间 碎片 的 轨道 特征 


四 理想 的 地 球 同步 轨道 是 圆 形 的 ， 具 有 轨道 倾角 ;= 0” ， 轨 道 平 运动 周期 为 1 d 等 特征 ， 这 样 的 轨 
道 也 称 为 静止 轨道 。 地 球 同步 轨道 高 度 为 353 785 km， 理 论 半 长 径 为 42 164 km。 事 实 上 ， 由 于 多 种 摄 
动因 素 的 影响 ， 地 球 同步 轨道 上 的 空间 碎片 存在 东西 向 和 南北 向 的 漂移 。 地 球 同步 轨道 的 主要 摄 动因 
素 为 地 球 的 非 球形 引力 摄 动 、 太 阳 和 月 亮 的 引力 摄 动 以 及 太阳 的 辐射 压力 摄 动 (5 。 太 阳 、 月 亮 的 引 
力 摄 动 造成 地 球 同步 轨道 目标 的 轨道 角 动量 有 53 a 左右 周期 的 进 动 ， 进 动 导致 轨道 倾角 存在 幅度 约 
为 15" 的 周期 变化 "1 ， 如 图 3。 地 球 同步 轨道 倾角 在 26. 5 a 后 会 由 0° 增 加 到 15” ， 再 过 26. 5 a 减 小 为 
0”。 太 阳 辐 射 压力 导致 一 个 小 的 偏心 率 ， 从 而 引起 径 向 半径 的 变化 ， 在 极端 情况 下 ， 这 个 变化 可 以 
达到 +75 km。 因 此 ， 定 义 地 球 同步 轨道 环 带 为 一 个 半径 为 和 2 164 km、 厚 度 为 150 km、 纬 度 + 15 的 区 
域 "( 如 图 4)。 


2 搜索 策略 分 析 与 计算 结 


对 于 地 球 同 步 轨 道上 的 空间 碎片 ， 即 使 以 理想 的 静止 轨道 为 对 象 ， 地 球 同步 轨道 环 带 区 域 里 的 空 
间 碎 片 实际 上 相对 地 面 测 站 是 运动 的 。 一 般 的 空间 碎片 观测 ， 每 个 目标 均 有 对 应 的 轨道 根 数 形式 的 轨 


@ http://www.scholarpedia.org/ article/ Space_debris 
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道 预 报 数据 。 设 备 根据 轨道 预报 数据 以 及 观测 站 svernal 


lequinox— 


的 位 置 ， 生 成 实时 和 望远镜 匹配 的 引导 数据 ， 引 
导 和 望远镜 指向 预报 位 置 ， 当 望远镜 到 达 预 报 位 置 
后 ， 由 于 导 星 镜 视 场 比较 大 ， 在 视 场 内 能 看 见 许 
多 星 ， 这 时 只 要 空间 目标 是 可 见 的 ,一般 可 以 依 
靠 观 测 人 员 的 经 验 对 目标 进行 人 眼 识 别 , 或 者 有 
一 些 识别 算法 ' 引 。 因 为 在 使 用 预报 数据 对 目标 
进行 跟踪 时 ， 目 标 在 视 场 内 相对 静止 ， 而 其 他 星 
是 运动 的 ， 以 此 可 以 对 目标 进行 判断 ， 随 后 选 定 
目标 ， 对 其 捕获 跟踪 “i( 即 凝 视 目标 )， 在 测 光 
视 场 内 对 其 进行 光度 测量 等 更 细致 的 观测 。 但 是 
当 目 标 非 常 瞳 弱 时 ， 在 导 星 镜 视 场 内 无 法 直接 看 
到 目标 ， 此 时 就 无 法 选 定 目 标 并 进行 跟踪 以 及 更 
细致 的 光度 测量 。 地 球 同步 轨道 空间 碎片 往往 属 
于 十 分 暗 弱 的 对 象 ， 尤 其 尺寸 相对 较 小 的 空间 碎 
片 ， 这 时 难以 直接 从 导 星 镜 视 场 里 找到 目标 ， 所 
以 需要 对 导 星 镜 视 场 进 行 搜索 。 Fig.3 Distributional relationship between orbital inclination 
由 于 预报 误差 ， 目 标 可 能 在 偏离 视 场 中 心 的 Et 

任何 一 个 位 置 。 大 假设 预报 数据 的 方位 角 精 度 
(标准 偏差 ) 为 o_A4， 高度 角 精度 为 o_EA， 并 假 
设 预 报 数据 与 目标 真实 位 置 的 误差 大 小 服从 高 斯 30*o2 100 km) 
分 布 。 从 而 真实 目标 偏离 预报 位 置 的 距离 值 也 为 02 150 km) 
高 斯 分 布 。 为 了 更 加 直观 地 分 析 搜 索 过 程 ， 将 预 
报 误差 角度 转换 到 目标 所 在 位 置 的 距离 值 (同样 
将 其 他 角度 值 转换 为 距离 值 ) : 

oI=Ro_A; 0o,= Ro_E, (1) 


其 中 , R 为 目标 与 测 站 的 距离 值 ，o,、o 分 别 
为 两 个 方向 上 偏离 距离 的 精度 经 过 上 式 转化 ， 相 
当 于 以 预报 位 置 为 中 心 ， 目 标 真实 位 置 离 中 心 位 
置 的 距离 ， 也 服从 高 斯 分 布 。 

以 预报 位 置 为 中 心 原点 ， 忽 略 径 向 ( 即 目 标 噬 离 ) 误 差 ， 建 立 一 个 目标 真实 位 置 偏离 值 的 二 维 高 
斯 分 布 1 ; 


图 3 ”地球 同 步 轨道 目标 轨道 倾角 (i) 与 升 交点 
杰 经 ( RAAN) 的 分 布 关系 


图 4 ”地球 同 步 轨道 环 带 区 域 
Fig.4 Geosynchronous orbital belt 
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其 中 ,原点 为 预报 位 置 ;，(x, y) 为 目标 的 真实 位 置 ; p(x, y) 为 目标 在 该 位 置 出 现 的 概率 ; o,、o, 为 
两 个 标准 差 ; pj,、j 为 两 个 均值 ; p 为 两 个 随机 变量 之 间 的 相关 系数 。 现 取 原 点 (预报 位 置 ) 为 中 心 ， 
人 i、Kt 则 均 为 0， 并 假设 两 个 方向 上 的 预报 误差 相互 独立 ， 因 此 相关 系数 为 0。 而 对 于 空间 碎片 ， 
大 部 分 目标 尤其 是 地 球 同 步 轨 道 碎 片 ， 由 于 缺乏 足够 的 观测 数据 ， 从 而 预报 精度 比较 低 ， 此 时 取 
ol=3000m，os=3000m， 计 算 目 标 真 实 位 置 在 距离 中 心 原点 20 km 内 出 现 的 概率 值 如 图 5、6。 

从 图 5、6 可 以 看 出 ， 越 靠近 预报 中 心 位 置 ， 目 标 出 现 的 概率 越 高 ， 越 远离 预报 中 心 ， 目 标 出 现 
的 概率 越 低 。 而 没有 颜色 的 区 域 并 不 代表 目标 不 会 出 现 ， 只 是 选取 了 圆心 为 零点 ， 以 1 km 为 半径 的 
区 域 进 行 了 高 斯 分 布 概率 计算 (该 区 域 全 概率 和 非常 接近 1) ， 圆 形 以 外 区 域 没 有 计算 ， 目 标 在 此 区 域 
外 出 现 的 概率 非常 低 。 
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图 5 目标 位 置 概率 分 布 ( 三 维 图 ) 图 6 目标 位 置 概率 分 布 (二 维 图 ) 
Fig.5 Target position probability distribution (3-D) Fig.6 Target position probability distribution (2-D) 


在 确定 目标 的 位 置 范围 以 及 分 布 概率 后 ， 为 了 成 功 地 搜索 到 目标 ， 需 要 在 该 范围 内 对 目标 进行 搜 
索 ， 一 般 采 取 的 搜索 方案 主要 是 同心 圆 搜索 (如 图 7) 和 把 楼 梯 搜 索 ( 如 图 8) 。 
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图 7 同心 圆 搜索 图 8 扑 楼 梯 搜 索 
Fig.7 Searching as concentric circles Fig.8 Searching as “climbing up” 


在 搜索 过 程 中 ， 考 虑 使 用 测 光 小 视 场 在 图 6 的 范围 区 域内 进行 搜索 ， 考 虑 测 光 视 场 的 大 小 ， 并 且 
通过 (1) 式 的 计算 ， 取 值 视 场 的 覆盖 范围 是 半径 为 1000m 的 区 域 。 从 而 为 了 找到 目标 ， 需 要 用 该 较 
小 的 视 场 对 目标 区 域 范围 进行 细致 的 搜索 ， 在 搜索 的 同时 适当 加 长 曝光 时 间 。 

图 7 即 是 同心 圆 搜索 ， 该 方案 从 搜索 区 域 中 心 位 置 开 始 ， 以 同心 圆 的 轨迹 一 次 次 向 区 域 边缘 搜 
索 ， 图 中 每 一 个 点 代表 一 次 搜索 ， 一 个 点 即 是 该 次 搜索 的 圆心 位 置 ， 覆盖 半 径 为 1000 m， 图 中 红线 
代表 搜索 轨迹 ， 顺 序 是 逆 时 针 方 向 。 在 计算 过 程 中 ， 已 经 保证 在 每 次 搜索 区 域 不 重 倒 。 对 于 第 n 次 搜 
索 ， 搜 索 覆 盖 范 围 可 以 表示 为 


(x -0) + (yy,) <R, (3) 
其 中 ，(%, ,yy ) 为 该 次 搜索 的 圆心 坐标 ; R 为 搜索 半径 ( 取 值 1000 m) 。 
搜索 成 功 的 概率 P 就 是 该 圆 形 区 域 蓝 盖 真实 目标 的 概率 ,结合 图 5 可 知 ， 通 过 在 圆 形 区 域 上 对 
(2) 式 表示 的 p(x, y) 进行 曲面 积分 ， 得 到 该 次 搜索 获 盖 目标 的 概率 值 : 
P=fp Cx, ydrdy , (4) 
积分 区 域 : (x -x,) ”+(y-y,) <R. 
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对 于 图 7 中 的 搜索 ， 从 圆心 开始 ， 沿 着 红线 ， 以 首 时 针 方 向 进行 同心 圆 搜 索 ， 对 每 次 搜索 进行 积 2 
计算 ， 得 到 单 步 搜索 成 功 的 概率 值 ， 可 以 得 到 如 图 9 的 结果 。 从 图 9 可 以 看 到 ， 每 步 的 搜索 概率 呈 阶 梯 
状 递减 ， 这 是 因为 圆心 在 同一 圈 上 的 搜索 区 域 对 应 的 概率 是 一 样 的 ， 随 着 同心 圆 搜 索 的 进行 ， 同 心 圆 半 
径 增 大 ， 处 在 同一 圆 上 的 搜索 区 域 的 概率 降低 。 在 图 9 中 ， 大 约 在 第 100 步 搜 索 时 ， 单 步 搜索 概率 就 
非常 接近 0， 后 面 的 搜索 概率 都 非常 低 ， 第 100 步 的 位 置 距离 中 心 原点 12 000 m， 即 说 明 12 000 m 半 
径 范 围 以 外 的 概率 已 经 非常 低 。 将 每 单 次 搜索 成 功 的 概率 累积 ， 就 可 以 得 到 图 10 中 虚线 代表 的 累积 
搜索 成 功 概率 ， 虚 线 表示 随 着 搜索 步 数 的 增加 ， 成 功 搜索 到 目标 的 概率 也 在 增加 。 结 合 图 9 与 图 10 
可 以 看 出 ， 大 约 在 第 100 步 搜 索 后 ， 单 步 搜索 概率 非常 接近 于 0， 累积 概率 增幅 几乎 可 以 忽略 。 
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Ql 让 楼 梯 搜 索 如 图 8， 以 红色 圆圈 标 出 起 始点 ， 
二 红色 线 代表 搜索 路 径 ， 蓝 色 点 代表 搜索 区 域 ， 同 
前 面 一 样 ， 图 中 每 个 点 是 一 个 圆心 位 置 ， 代 表 一 证 
个 半径 为 1000m 的 圆 形 区 域 。 以 同样 的 方式 ， 在 
每 个 点 代表 的 搜索 区 域 进行 积分 运算 ， 按 顺序 计 yw 004 
[= 算 每 个 搜索 区 域 的 成 功 概率 如 图 11。 是 003 
OQ 图 11 每 次 搜索 成 功 的 概率 随 着 搜索 步 数 变化 ， 守 
结合 图 8 中 的 搜索 路 径 可 知 ， 搜 索 区 域 与 中 心 原点 。 3 0.02 
的 距离 值 从 小 到 大 反复 变化 ， 从 而 造成 了 图 11 中 本 
概率 值 的 变化 。 同 样 ， 将 每 步 的 搜索 概率 一 步 步 累 
上 ， 从 而 得 到 图 10 中 的 红色 实 线 ， 该 线 表示 随 着 搜 0 
索 步 数 的 增加 ， 成 功 搜索 到 目标 的 概率 也 在 增加 。 0 
从 图 10 可 以 看 出 ， 两 种 方案 虽然 搜索 了 同样 
的 步 数 ， 但 两 方案 的 累积 搜索 概率 差别 较 大 ， 两 人 


Fig. 11 Successful search probability of each step 


种 方案 同样 一 共 搜索 了 271 步 。 蓝 色 虚 线 表示 的 
同心 圆 方案 在 大 部 分 搜索 步 数 时 保持 着 更 高 的 搜索 成 功 概 率 ， 并 且 此 方法 可 以 很 快 达到 概率 峰值 ， 也 
就 是 搜索 效率 比较 高 ， 该 方案 在 第 94 步 累积 概率 值 为 0.734 0， 随 后 在 129 步 达到 概率 峰值 0. 734 8。 
拒 楼 梯 搜 索 方案 初始 的 步 数 概率 很 低 ， 随 后 很 快 升 高 ， 最 终 在 第 178 步 ， 疏 楼 梯 搜 索 成 功 概 率 刚 
好 等 于 同心 圆 搜索 概率 (0. 734 8) ， 随 后 该 方法 在 第 219 步 累 积 概率 达到 0. 784 4， 最 终 候 楼 梯 方 法 在 
第 250 步 达 到 概率 峰值 0. 785 4。 
由 以 上 分 析 可 知 ， 同 心 圆 搜 索 方 法 搜索 效率 更 高 ， 可 以 在 较 少 的 步 数 之 内 达到 较 高 的 搜索 成 功 概率 。 


3 期 汤 侍 峰 等 : CEO 暗 弱 空间 碎片 目标 搜索 策略 分 析 309 


3 结论 与 展望 


不 考虑 预报 精度 在 方位 角 与 高 度 角 差 异 的 情况 下 ， 采 取 同 心 圆 搜索 方案 比较 高 效 ， 由 文中 的 计算 
可 知 ， 该 方案 在 第 94 步 达 到 较 高 的 累积 概率 ， 而 往 后 的 搜索 步 数 效 率 非常 低 ， 可 以 考虑 开始 新 的 一 
轮 搜索 。 实 际 搜 索 中 ， 在 大 致知 道 预报 误差 的 情况 下 ， 可 以 根据 文中 方法 估算 最 适合 搜索 的 区 域 ， 从 
而 提高 搜索 效率 。 而 当 预 报 误差 精度 不 一 致 时 ,通过 合理 地 选取 起 始 搜 索 点 ， 应 该 可 以 提高 候 楼 梯 搜 
索 的 效率 ， 甚 至 可 以 使 得 此 方法 更 加 优秀 。 
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Abstract: The geosynchronous region is one of the most valuable space environment for application which 
makes it of great significance to monitor the space debris by optical observation. On the basis of the special 
dynamical characteristics of geosynchronous orbits, the motion features of CEO space debris were introduced in 
this paper. Concerning the prediction error of the space target and size of field of view, the success 
probabilities of “ concentric circles search”and“climb stairs search” were calculated. Based on the calculated 
results, the strengths and weaknesses of each search strategy were analyzed. According to the analysis, the 
“ concentric circles search” is more efficient in most cases. 
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